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ABSTRACT 

We present spectroscopy and laser guide star adaptive optics (LGSAO) images of 
the doubly imaged lensed quasar SDSS J1206+4332. We revise the deflector redshift 
proposed previously to zg = 0.745, and measure for the first time its velocity dispersion 
a = (290 + 30) km/s. The LGSAO data show the lensed quasar host galaxy stretching 
over the astroid caustic thus forming an extra pair of merging images, which was 
previously thought to be an unrelated galaxy in seeing limited data. Owing to the 
peculiar geometry, the lens acts as a natural coronagraph on the broad-line region 
of the quasar so that only narrow [O III] emission is found in the fold arc. We use 
the data to reconstruct the source structure and deflector potential, including nearby 
perturbers. We reconstruct the point-spread function (PSF) from the quasar images 
themselves, since no additional point source is present in the field of view. From 
gravitational lensing and stellar dynamics, we find the slope of the total mass density 
profile to be 7 ' = — log p/log r = 1.93 ± 0.09. We discuss the potential of SDSS 
J1206+4332 for measuring time delay distance (and thus Ho and other cosmological 
parameters), or as a standard ruler, in combination with the time delay published by 
the COSMO GRAIL collaboration. We conclude that this system is very promising 
for cosmography. However, in order to achieve competitive precision and accuracy, 
an independent characterization of the PSF is needed. Spatially resolved kinematics 
of the deflector would reduce the uncertainties further. Both are within the reach of 
current observational facilities. 

Key words: gravitational lensing: strong - quasars: emission lines - galaxies: kine¬ 
matics and dynamics - cosmology: distance scale - methods: observational - methods: 
statistical 


1 INTRODUCTION 

In the past two decades, the dream of Refsdal (1964) to 
use gravitational time delays as a tool for cosmography has 
become a reality. Several observational and modeling break¬ 
throughs now make it possible to measure time delay dis¬ 
tances with 5-6% precision and accuracy from a single galaxy 
scale gravitational lens system (e.g. Suyu et al. 2010), and 
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the angular diameter distance to the deflector (Parahcz & 
Hjorth 2009; Jee et al. 2014). In turn, these absolute dis¬ 
tance measurements enable the determination of the Hub¬ 
ble Constant (e.g., Schechter et al. 1997; Treu Sz Koopmans 
2002; Wucknitz et al. 2004; Jackson 2007; Oguri 2007), and 
other cosmological parameters like flatness and equation of 
state of dark energy when combined with other probes (Suyu 
et al. 2013, 2014). The necessary observational ingredients 
include determination of time delays with few percent pre¬ 
cision (Burud et al. 2002; Fassnacht et al. 2002; Eulaers et 
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al. 2013; Tewes et al. 2013), imaging and astrometry at 10- 
100 mas resolution (Fassnacht et al. 2002; Koopmans et al. 
2003; Chantry et al. 2010), stellar velocity dispersion of the 
deflector (Treu & Koopmans 2002), and imaging and spec¬ 
troscopy to characterize the line of sight and environment 
(Keeton &; Zabludoff 2004; Fassnacht et al. 2011; Greene et 
ah 2013; Collett et ah 2013). 

The application of lensed quasars for cosmography has 
been limited by the small number of systems for which the 
full complement of data is available. We present here results 
from on ongoing program, aimed at increasing the statistical 
power of the method, by obtainining all the required imaging 
and spectroscopic data for gravitational lensed quasars with 
measured time delays. 

In this paper, we focus on a very unusual and remark¬ 
able system, the doubly lensed quasar SDSS J1206+4332. 
It was discovered by Oguri et al. (2005), who suggested a 
redshift of z = 0.748 for the lens galaxy based on Mg II 
absorption on the quasar spectra, and has been monitored 
by Eulaers et al. (2013) to obtain a time-delay At = 111 ±3 
days between the two quasar images. In those observations, 
an additional galaxy ‘G3’ was detected along the line of 
sight NE of the system, which would complicate the lens 
modelling. Our laser guide star adaptive optics (LGSAO, 
Wizinowich et al. 2006) data show that ‘G3’ is in fact two 
merging images of lensed quasar host galaxy itself. In prac¬ 
tice, while the quasar is just outside the astroid caustic and 
thus only doubly imaged, the host galaxy cross the caus¬ 
tics and produces four images in a fold configuration. This 
configuration makes this system very interesting for cosmog¬ 
raphy, since it is relatively easy to monitor as a double, and 
yet it has a larger number of imaging constraints available 
for reconstructing the lensing potential than a typical pure 
double. 

We complement our LGSAO data with long-slit spec¬ 
troscopy of the deflector, whence we can revisit the red- 
shift and provide the first measurement of the velocity 
dispersion. Spectroscopy and high-resolution imaging en¬ 
able the reconstruction of source structure and inference 
of the gravitational lens potential. Together with the pub¬ 
lished time-delay, they enable us to assess the role of SDSS 
J1206+4332 as a probe of the Hubble constant, or alterna¬ 
tively as a standard ruler. 

This paper is structured as follows. In Section 2 we de¬ 
scribe the imaging data, complemented with the spectro¬ 
scopic follow-up in Section 3. In Section 4, we illustrate 
the reconstruction of deflector, source and PSF. We add 
information form stellar kinematics in Section 5, obtain¬ 
ing an independent characterization of the deflector total 
mass density profile. We then discuss the value of SDSS 
J1206+4332 for cosmography, in terms of inference on Ho 
and deviations from the flat ACDM paradigm. We conclude 
in Section 7. 


2 KECK ADAPTIVE-OPTICS IMAGING 

We imaged SDSS J1206+4332 in the K' band at 2.2/xm with 
the NIRC2 camera on Keck II on 2014 January 8th, using 
the LGSAO system. A total exposure time of 4860 s was 
obtained with 27 exposures of 180 s each. Data were reduced 
and coadded with pipelines developed by Auger et al. (2008). 
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Figure 1 . Imaging of SDSS J1206+4332 with LGSAO-NIRC2; 
North is up and East is left, the bar marks 1". The two quasar 
images, the deflector galaxy (‘Gl’ in Oguri et al. 2005; Eulaers et 
al. 2013) and the lensed host galaxy are all clearly visible. The 
lensed host traces a typical fold configuration, previously dubbed 
£ G3’ and mistaken in the past for an additional galaxy along the 
line of sight. Three faint objects are visible to the NW, at w 4.4 X/ 
from the deflector, dubbed £ G2’ in previous observations where 
they were blended into a single, broad and faint object. The green 
horizontal lines mark the extent of the DEIMOS slit used for 
spectroscopy. 


High resolution is obtained using the pixel scale of 9.94 x 
10 -3 arcseconds of the Narrow Camera of NIRC2. 

The coadded data are shown in Figure 1, rebinned in 
three-by-three pixels to highlight low surface-brightness fea¬ 
tures. Here and through the rest of this paper, North is up 
and East is left. With this resolution, the two quasar im¬ 
ages are clearly visible and distinguished from the deflector 
galaxy in the middle. The (lensed) host galaxy is visible 
around the quasar images and also in the fold arc at NE, 
which was previously mistaken for a distinct blue galaxy 
(dubbed ‘G3’; Oguri et al. 2005; Eulaers et al. 2013). 

We fit the deflector light with a flattened Sersic pro¬ 
file, within a circular patch avoiding contamination from 
the Einstein ring light. To this aim, the model profile must 
be convolved with the point-spread function (PSF). Given 
the small field of view, there is no independent point source 
at our disposal for this purpose, so we must infer the PSF 
directly from the system itself. This is described in the sub¬ 
section 2.1. The resulting parameters of the deflector light 
are given in Table 1. We will use the best-fitting light-profile 
in the following Sections to model the lensing and dynamics. 
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Figure 2. Cutouts around the two point-sources (top line), their 
scaled difference (bottom left) and first guess fo the PSF (bottom 
right) as decribed in Section 2.1. 


b/a 

p.a. (N 

of W) 

^ Reff 

Ara 

Adec 

0.89 

51.24 

deg 

2.66 0.60" 

0.252" 

2.583" 


Table 1. Parameters of the best-fitting deflector light profile. 
From left to right: axis ratio, position angle (counterclockwise, N 
of W); Sersic index; effective radius; relative astrometry from the 
southern-most quasar image, peak to peak, with r.a. (resp. dec.) 
increasing to the East (resp. North). 


2.1 PSF estimation 

For the deflector light subtraction, a first guess of the PSF 
is sufficient. We obtain it by isolating small cutouts centred 
on the two quasar images, subtracting to each of them the 
average flux around the border, normalizing their fluxes to 
unity, then coadding the two cutouts. 

To recover a PSF that is common to both images and 
avoid propagating noise, we use a crude version of regu¬ 
larization. For each pixel of the target PSF coadd we con¬ 
sider the two corresponding pixel-values p ±, P 2 in the point- 
source cutouts: when these are within 0.1 of the average 
av p — (Pi +P2)/2, we register av p in the PSF coadd, oth¬ 
erwise we register the minimum min(pi,p 2 )- The result is 
shown in Figure 2. 

The reconstruction relies on the fact that the lensed 
host is stretched tangentially by gravitational lensing, so 
that in first approximation it contributes a uniform back¬ 
ground flux to the point-source, which should dominate over 
the PSF light at large enough distances from the PSF core. 
A variant of this procedure will be combined iteratively with 
the lens model in Section 4. 


3 KECK SPECTROSCOPY 

Long-slit spectroscopic measurements of SDSS 
J1206+4332 were taken with the instrument DEIMOS 
(Faber et al. 2003) on Keck 2, on May 16 2015. The l”-wide 
slit was centred on the dehector galaxy and aligned with 



Observed wavelength (A) 


Figure 3. DEIMOS slit spectra. Top: Extracted ID spectrum 
of the quasar host galaxy. Bottom: Extracted ID spectrum of the 
quasar (northern point source) and the dehector galaxy. The red 
line is the best-ht spectrum obtained from the velocity dispersion 
fitting. The green line is the best-ht polynomial used to model 
continuum emission from the quasar. The vertical gray bands are 
regions of the spectrum masked out from the ht. 


the long axis in the E-W direction (see Fig. 1). 

We used the 600ZD grism, covering the wavelength range 
4600A - 9200A , with a spectral resolution of about 
160 km s -1 FWHM. The total integration time was 1 hour. 
We detect a signal from the dehector galaxy, the quasar, 
and the quasar host. The lens galaxy and quasar host are 
spatially resolved, while the quasar itself, not covered by 
the slit, shows up as a contamination to the lens galaxy 
spectrum. The ID spectra of the lens galaxy and quasar 
host are plotted in Figure 3. 

The quasar component is visible with its continuum 
emission in the blue side of the lens galaxy spectrum and 
with broad emission lines. The host galaxy spectrum is very 
faint. The only clear line detected is [C III] 1908A, also seen 
in the quasar component. However, while the quasar line is 
broad (FWHM ~ 6000 km s -1 ), the same line is much nar¬ 
rower in the spectrum of the host. Thus, we infer that the 
quadruply imaged part of the host galaxy does not include 
the compact broad line region, but only the significantly 
more extended narrow line region (NLR). In other words, 
the lens is acting as a natural “coronagraph”, blocking out 
the light of the accretion disk and the BLR in two of the 
images. From a ht to the narrow line wavelength, we obtain 
a source redshift = 1.789, confirming the one measured 
by Oguri et al. (2005). At lower signal-to-noise ratio, there is 
an absorption feature compatible with Fell at 2344A at the 
host redshift z s , which is seen also in the quasar spectrum. 

The lens galaxy spectrum has some prominent stellar 
absorption features: Ca K,H at 3934A, 3967A, the G-band 
absorption complex around 4300A, and Mgb at 5175A, at 
a redshift Zd — 0.745. In addition to these stellar absorp¬ 
tion lines we detect nebular absorption in the Mg II doublet 
2795A - 2803A at z — 0.748. These lines were also detected 
by Oguri et al. (2005), who used them to estimate the lens 
redshift. There is a rest-frame velocity difference of about 
516 km s -1 between stellar and Mg II absorption. We note 
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that this difference in redshift would introduce a small sys¬ 
tematic bias in cosmological measurements (e.g. 0.15% in 
angular diameter distance to the deflector) based on this 
system if not corrected. 

From the data at hand we can also extract a velocity 
dispersion measurement. We do this by fitting stellar tem¬ 
plates to the observed spectrum, using an adaptation of the 
velocity dispersion fitting code by van der Marel et al. (1994) 
described by Suyu et al. (2010). We also fit for continuum 
emission from the quasar with a 7-order polynomial. We use 
G and F star templates from the fndo-US library. The ID 
spectrum of the lens is extracted from a l”xl.9” rectan¬ 
gular aperture. We fit the rest-frame wavelength region be¬ 
tween 3800A and 5250A, where the most important stellar 
aborption features are found, masking out the region con¬ 
taminated by broad Mgll emission from the quasar as well 
as the deep atmospheric feature around 7600A. The deepest 
stellar absorption lines, Ca K,H, fall in correspondance with 
the telluric absorption feature at 6900A. We correct for this 
atmospheric absorption using a telluric standard star ob¬ 
servation as a reference and keep Ca K,H in the region of 
the spectrum used for the fit. The best-fit spectrum is over- 
ploted in Figure 3. The median signal-to-noise ratio in the 
region used for the fit is 2lA _1 , though if we consider only 
the signal coming from the lens galaxy, i.e. once contami¬ 
nation from the quasar continuum emission is removed, it 
decreases to 13A -1 (rest-frame). 

We measure a velocity dispersion of a = 290 db 
30 km s- 1 . The uncertainty reflects both statistical and sys¬ 
tematic errors. We estimated the latter by repeating the fit 
over different wavelength ranges and by changing the order 
of the polynomial continuum. 


4 LENS MODEL 

We use the modelling code GLEE (Suyu & Halkola 2010; 
Suyu et al. 2012 ) to infer the lensing potential and recon¬ 
struct the extended source structure. As constraints, we use 
the arc surface-brightness distribution and the point source 
positions. The arc is cleaned from possible contamination by 
the lens light, by subtracting the best-fitting profile found in 
Section 2. Due to uncertainties in the PSF reconstruction, 
we artificially increase the pixel uncertainties around the 
point-source images, so that the lens model is constrained 
mainly by the NE fold and by the two quasar image posi¬ 
tions, to which we assign a positional uncertainty of 20 mas. 
For computational simplicity, we consider just the pixels in 
an annular mask around the Einstein radius. 

Inference on the lens mass profile relies primarily on 
the differential distortion and radial magnification of the 
same source patch on different sides of the lens. Then, two 
conditions must be met: (i) high resolution imaging, with 
good seeing conditions and small pixel size; (ii) a robust 
characterization of the PSF. While the former is satisfied by 
our Keck images, the latter is not directly available and the 
PSF must be reconstructed from the images. 

The full model would then encompass both the lens 
parameters and the PSF reconstruction. We proceed itera¬ 
tively, alternating PSF reconstruction and lens model fitting. 
From a run of the lens model, we get the predicted surface- 
brightness profile of the lensed host, which we can subtract 


near the quasar images. Similarly to the procedure in Sec¬ 
tion 2 . 1 , we select two cutouts around the (host-subtracted) 
point-source images and combine their common pixel-values 
into a new PSF, which is then used in a new run of the lens 
model. Figure 4 shows the inferred PSF in subsequent iter¬ 
ations. After three iterations, the correction step starts to 
overfit noise, while the parameters of the lens model are not 
changing appreciably. 


4.1 Model Components 


We adopt a power-law profile for the total (i.e. luminous 
plus dark matter) mass density of the main deflector, with 
convergence 


n{x,y) 


3-7' / yr 2 + g-y 

1 + q y Re 


i-Y 


a) 


along the principal axes of the deflector. With this conven¬ 
tion, the Einstein radius Re is such that in the circular limit 
(q = 1 ) it encloses a mean convergence of 1 . In what follows, 
the model Einstein radii will be in angular units and denoted 
by 0 E - 

The main deflector has a three nearby, faint perturbers 
to the NW, at a distance of approximately 4.4”, as shown 
in Figure 1. tn previous observations (Oguri et al. 2005; Eu- 
laers et al. 2013), they were all blended into a single, broad 
and faint clump dubbed C G 2 ’. In order to model contribu¬ 
tions from the immediate environment, we then include a 
Spherical Isothermal Sphere (SIS, 7 ' = 2 , q — 1 ) centred on 
G 2 , with a centre to be adjusted, plus external shear. The 
role of shear is to correct for: flattening and/or imperfect 
centering of the group G 2 ; lesser contributions from other 
galaxies in the vicinity of the lens; mismatch between the 
adopted and real flattening of the main deflector. 

The model then has the following parameters: the 
power-law exponent 7 ' of the density profile of the deflector, 
its axis ratio q and position angle (pi (N of W); the shear 
amplitude and position angle ip s (N of W); the nearest 
perturber (G2) Einstein radius 6e, p and its position; the 
source-position, constrained by the quasar image positions 
within a 20 mas range, and the point-source fluxes /i,/ 2 - 
The point-source fluxes are computed but not used for the 
lens model, as they can be affected by source variability, 
microlensing, and milli-lensing by possible satellites of the 
main deflector. One final parameter would the Einstein ra¬ 
dius 6e,i of the main deflector, which we keep fixed because 
it can be determined with 1 % precision as shown below. 


4.2 Initialization and Choice of Parameters 

initialization is performed via conjugate-point analysis, i.e. 
identifying sets of points that should map to the same loca¬ 
tion in the source-plane and fitting for their predicted po¬ 
sitions in the image-plane. Given the low number of con¬ 
straints, we can examine just combinations of separate mod¬ 
els, whose parameters will be cross-checked for consistency. 
The first initialized model has an isolated Singular fsother- 
mal Ellipsoid (StE, 7 ' = 2, q < 1), with free Einstein radius, 
axis ratio and position angle. The second model has a StS 
plus external shear, the third has a StS for the main deflector 
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model 

0e,i 

(arcsec) 

q 

Vl 

[rad] 

7s 

Vs 

[rad] 

#E,p 

(arcsec) 

Ara 

(arcsec E) 

Adec 
(arcsec N) 

SIE 

1.61 

0.807 

0.65 

— 

— 

— 

— 

— 

SIS+XS 

1.582 

— 

— 

0.05 

1.95 

— 

— 

— 

SIE+XS 

1.582* 

0.814 

0.87 

0.04 

1.60 

— 

— 

— 

SIS+SIS 

1.571 

— 

— 

— 

— 

0.30 

-2.506 

3.485 

SIE+SIS 

1.571* 

0.847 

0.485 

— 

— 

0.25 

-2.198 

3.852 

all 

1.575* 

0.8* 

0.75 

0.034 

1.12 

0.25 

-2.352* 

3.669* 

all 

1.575* 

0.86* 

0.82 

0.025 

1.56 

0.18 

-2.352* 

3.669* 

all 

1.575* 

0.9* 

0.91 

0.032 

1.89 

0.17 

-2.352* 

3.669* 

all 

1.575* 

0.99* 

1.10 

0.05 

2.2 

0.16 

-2.352* 

3.669* 


Table 2. Model parameters obtained from conjugate-point analysis, using different models. Parameters marked by an asterisk are kept 
fixed during the optimization. These are then used to initialize a model with main lens, external perturber G2 and external shear, which 
is required to fit the quasar image positions and arc surface-brightness profile. The Einstein radius of the main lens is kept fixed to its 
SIS value when generalizing to SIE. The second part of the table shows results for a model with SIE, external shear and perturber, with 
some paramters fixed. 



Figure 4. Reconstructed lens model. Top: iterative corrections of the reconstructed PSF. A 23-by-23 pixel cutout was chosen for the 
PSF; since the third iteration, the small-scale PSF corrections start to overfit noise in the data. Bottom: deflector-light-subtracted data 
(left), model of the arc (middle) and residuals over noise (right). Secondary lobes in the PSF are evident in the deflector-light-subtracted 
image and are masked out when fitting the extended structure. Secondary ‘wings’ in the reconstructed profile around the quasar images 
(near the mask border) are given by residuals in masking the PSF lobes out and choice of error-map. Their appearance is also degenerate 
with the point-source fluxes, which can create ‘holes’ in the reconstructed arc. The residuals on the fold arc are completely compatible 
with noise. 


and one for G2. Since the SIS models yield consistent mea¬ 
surements of Einstein radius, we then generalize the main 
deflector to a SIE, keeping its Re to the SIS value, and op¬ 
timize a model with external shear and one with a SIS for 
the external perturber G2. 

The results are shown in Table 2. The position of G2 is 
consistent within « 0.4 // . Its inferred inferred values of Oe 
are consistent within 0.05”. The difference in #e,i between 
a model with external shear and a model with external SIS 
perturber is ~ 0.01 // . 


The interplay between different parameters is explored 
by different choices of flattening of the main deflector, keep¬ 
ing the perturber position fixed to the average value between 
the SIS+SIS and SIE+SIS cases. Results from this experi¬ 
ment are shown in the second half of Table 2. When q — 0.9 
as derived from the photometry of the deflector galaxy in 
Sect. 2, the inferred deflector p.a. from lensing is exactly 
the same as the one found from the photometry. This find¬ 
ing suggests that q — 0.9 and ipi — 0.91 can be kept fixed 
in the lens model. Similarly, we keep the perturber position 
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Figure 5. Reconstructed lens model. Inference on the four main 
parameters of the lens model. We show the marginalized posterior 
distribution of the power-law slope 7 ', external shear amplitude 
7 s and angle (/? s , and the perturber (G2) Einstein radius #E,p i n 
arcseconds. Due to the very small pixel-size and lack of indepen¬ 
dent PSF, the model uncertainties are artificially small. 

fixed, since possible corrections will be given by the shear 
term. 

Unfortunately, with the data at our disposal, a complete 
lens model on the extended source is unstable and some of 
the parameters converge to pathological minima. For this 
reason, we are currently forced to keep some of the param¬ 
eters fixed to their fiducial values. 


4.3 Lensing Results 

Results from the model are shown in Figure 5, showing the 
inferred posterior distribution on the four main parameters. 
Their values are in general agreement with what is found 
at the initialization stage. The Einstein radius 0e, p of G2 is 
small but non negligible. Other perturbations to the lensing 
potential are summed up in a small shear, which corrobo¬ 
rates our working hypothesis. The density slope 7 ' = 2.04 
is in line with the general description of early-type galaxy 
lenses (Koopmans et al. 2009), but steeper than the median 
value ~ 1.95 of lens galaxies at this redshift (Sonnenfeld et 
al. 2013). 

The statistical uncertainties are driven by the small 
pixel size compared to the Einstein radius, and are therefore 
dominated by systematics in the choice of lens model, PSF 
correction scheme, and estimation of errors near the point- 
source images. Results on other systems with HST data, 
with a robustly characterized PSF and deeper observations 
on the Einstein ring, typically yield uncertainties around 
2 — 3% on 7 ' (Dye & Warren 2005; Suyu et al. 2013). In 
what follows, we will adopt an estimated 5% for the uncer¬ 
tainties on 7 '. 

Galaxies in the vicinity of the main deflector contribute 
to the deflections of light-rays. As a first approximation, 
their role can be subsumed in a uniform sheet of conver¬ 
gence K e around the main deflector. However, this can be 
reabsorbed by (unobservable) shifts and rescalings in the 
source-plane, so that K e cannot be determined from lensing 



Figure 6. Combined inference from lensing and dynamics. The 
model velocity dispersion is plotted against total density expo¬ 
nent 7 ' = —d log p/d log r, corresponding to the measured Ein¬ 
stein radius $e,i = 1.575". We also show the measured velocity 
dispersion with its uncertainty (horizontal lines) and the density 
exponent from the lens model, with a conservative uncertainty of 
5% (vertical lines). 


alone. This mass-sheet degeneracy (MSD, Falco et al. 1985) 
is a simple case of source-position transformation degenera¬ 
cies (Schneider &; Sluse 2013), with particular importance 
for time-delay cosmography. Bounds on K e can be put by 
studying the environment of the lens (Fassnacht et al. 2011; 
Greene et al. 2013). Oguri et al. (2005) have shown that the 
environment of SDSS J1206+4332 is comparable to a ran¬ 
dom field of view, considering galaxies within 1 arcmin from 
the lens down to an /—band magnitude of 24. Then, given 
the current data quality, we can neglect n e for our purposes. 
We will discuss this choice further in the following Sections. 


5 DYNAMICS 

Information on the kinematics of stars in the deflector yields 
and independent mass probe, thus breaking residual degen¬ 
eracies in the lens model (Treu & Koopmans 2002, 2004). 
This is especially true for the power-law profiles used here, 
for which simple relations hold between the velocity dis¬ 
persion and power-law exponent (Koopmans 2006). From 
gravitational lensing, a normalization condition holds for the 
projected mass of the lens within the Einstein radius: 

2 T^\ 

M p (R E ) = (1 - Ke)7Tj?|S cr = JL S ~ (1 - K e )R 2 E . ( 2 ) 

4G JJi J-sis 

In terms of the mass M(r), surface-brightness E(R) = 
£o Ti(R/Re) and 3D luminosity v(r) — i>(r/R e )^o/R e pro¬ 
files, the line-of-sight velocity dispersion a p (R) satisfies 

e (R)a 2 p (R) = 2 g f°° ( y r 2 _ + kp ( Ttfl )) dr 

Jr r 

(3) 

where the kernel kg depends on the orbital anisotropy pro¬ 
file /3(r) (Mamon & Lokas 2005; Agnello et al. 2014). In the 
power-law case, with p(r) oc r -7 , the lensing mass normal- 
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ization yields 


E (R/R e )<r 2 p (R) 


(1 — K e )c 2 


D s Re 

DisDi 


x 


(3-7 , )r( 7 V2) 
\/*T((7' - l)/2) 


X 



(V® 2 - (R/Re) 2 + kp(x, R/Rej) d* 


The model velocity dispersion can be averaged over the 
slit or aperture to be directly compared with its measured 
value, either via numerical integration or with the more com¬ 
pact formalism of Agnello et al. (2013), which is adequate 
here because most of the deflector light is enclosed in the 
DEIMOS slit. Then, given 0e = Re/Di and the distance- 
ratio D s /Di s , from the measured kinematics one can derive 
the exponent 7 ' and vice versa. In Figure 6 we show the pre¬ 
dicted velocity dispersion as a function of 7 ', together with 
the measured value and the density exponent inferred from 
the lens model. 

Figure 6 shows the predicted velocity dispersion a as 
a function of 7 '. The vertical lines delimit the 68 % confi¬ 
dence interval in 7 ' predicted by the lens model, while the 
horizontal lines show the measured a with its 68 % confi¬ 
dence interval. The data are not sufficient to discriminate 
between models with different anisotropy, since we just have 
an aperture-averaged value (with broad PSF) instead of a 
radial profile at our disposal. 

The slope obtained this way, 7 'm 1.93 ±0.09, is smaller 
than the value inferred from lensing alone. Still, when a con¬ 
servative uncertainty of 5% on the lensing 7 ' is adopted, the 
two results are compatible within 68 % confidence level. 


6 PROSPECTS FOR COSMOGRAPHY 

The time-variability of the source quasar provides an addi¬ 
tional observable, which is the time-delay between the two 
quasar images. The arrival time is given by 

t = ( 1 + z d)^^- [(x;-x s ) 2 -^(x;)] (5) 

in terms of the projected lens potential - 0 , deflection x/ -xg 
and cosmological distances. 

From the lens model, the predicted time-delay can be 
compared with its measured value to infer the time-delay 
distance D/± t = (1 ± Zd)DiD s /Dd s oc Hq 1 , or equivalently 
the cosmological parameters, primarily the Hubble constant. 
This is generally plagued by the MSD, as the source-position 
transformation that eliminates K e from the lens model pro¬ 
duces a biased measurement Ho /(.1 — K e ) of the Hubble 
constant, with everything else fixed - cosmological abun¬ 
dances and measured At. If the velocity dispersion and 
time-delay are modelled jointly, the common dependence on 
(1 — K e )D s /D ds cancels out, so that the angular-diameter 
distance Dd to the deflector is inferred directly, with possi¬ 
ble systematics from K e now erased (Jee et al. 2014). Alter¬ 
natively, in presence of a good prior on K e , dynamics can be 
used to obtain the exponent 7 ' to be used in the time-delay 
inference, since the model kinematics depend just weakly on 
cosmological parameters (eq. 4). 

Given the data presented here, with small nominal un¬ 
certainties on 7 ', one would be tempted to exploit them 
for time-delay cosmography. We caution against that, for 
the following reasons. First, the signal-to-noise ratio of faint 


features, such as the fold arc and perturbers in G 2 , is too 
low to enable a flexible model with robustly determined pa- 
(4) rameters. Second, the lack of an independent PSF model 
leaves degeneracies in the lens model unresolved, degrading 
the predictive power near the point-source positions. This 
is particularly true for two-quasar-image systems like SDSS 
J1206±4332; in four-quasar-image systems, there is more 
information from the four quasars to reconstruct the PSF 
directly from the lens system (Chen et al., in preparation). 
The confidence ranges of lensing parameters account just 
for the statistical uncertainties and are driven mainly by 
the small pixel-size, whereas systematic uncertainties (e.g. 
from PSF correction) dominate the error budget. 

Bearing the caveats in mind, it is interesting to provide 
some rough cosmographic estimates and discuss what can be 
earned from more reliable data. For concreteness, we forecast 
the cosmographic constraints that we expect to obtain with 

(1) HST imaging that have better characterized PSFs, and 

( 2 ) spatially resolved kinematics of the deflector galaxy that 

help break lens-modeling and mass-sheet degeneracies. To 
estimate the precision on the time-delay distance, we assume 
uncertainties on the time delays of 3% (Eulaers et al. 2013), 
modeling with HST imaging of 3% (e.g., Suyu et al. 2013), 
and external convergence K e of 4% (e.g., Greene et al. 2013, 
Collett et al. 2013). Adding these contributions in quadra¬ 
ture, we forecast an uncertainty of 5.8% on R>At • For refer¬ 
ence, we adopt a fiducial value for D/± t which corresponds 
to a fiducial flat ACDM model with fi m — 1 — — 0.3 and 

Ho = 70 kms - 1 Mpc -1 . Assuming a Gaussian probability 
distribution for D/± t — 5789 ± 336 Mpc, we show in Figure 7 
the cosmographic information from SDSS J1206±4332 as¬ 
suming Planck priors (Planck Collaboration et al. 2014) in 
either open ACDM cosmology (left panel) and flat reCDM 
cosmology (right panel). Such a distance measurement of 
DAt for SDSS J1206±4332 would help break significantly 
parameter degeneracies in the CMB data, and provide an 
independent test of systematics in cosmographic probes. 


7 SUMMARY 

We have presented laser guide star adaptive optics images 
and spectroscopic observations of the gravitational lens sys¬ 
tem SDSS J1206±4332. A precise time delay measurement is 
available for this doubly imaged quasar (Eulaers et al. 2013), 
making it a prime target for follow-up observations aimed at 
using it as tool for cosmography. The high resolution images 
allow us to recognize that the system is in a very unusual 
and favorable configuration for cosmography, previously not 
resolved in seeing based images. Furthermore, our deep Keck 
spectroscopy allows to correct the previous estimate of the 
redshift of the main deflector and measure its stellar velocity 
dispersion based on stellar absorption features. We then use 
the newly available information to construct a gravitational 
lens model of the system and discuss its potential for cos¬ 
mography. Our main results can be summarized as follows: 

• The quasar lies just outside of the astroid caustic in 
the source plane and is thus doubly imaged. However, the 
quasar host galaxy crosses the caustic and produces two ad¬ 
ditional merging images in a classic fold configuration - pre¬ 
viously mis-identified as an unrelated galaxy along the line 
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Figure 7. Predicted inference on cosmological parameters, in presence of robust model parameters (with well controlled systematics) 
from more reliable data than the ones available here. Dashed red lines are the Planck prior on cosmological parameters (Planck Collab¬ 
oration et al. 2014), the black lines are the marginalized posterior using SDSS J1206+4332 as a cosmographic probe, assuming state of 
the art follow-up data are available. Degeneracies between Ho and other cosmological parameters are appreciably reduced by time-delay 
lenses, which probe primarily Ho. Left: Open ACDM cosmology (free curvature parameter Qk)- Right: Flat CDM (fife = 0) with free 
equation of state parameter w. 


of sight in seeing limited images. The presence of two addi¬ 
tional images doubles the amount of information available 
to constrain the gravitational lens model. 

• The lensing geometry acts as a natural “coronagraph”, 
blocking out the light coming from the broad line region 
and the accretion disk in two of the images. Thus, those 
two images give a very clear view of the narrow line region 
surrounding a z = 1.79 quasar. This system represents an 
interesting target for follow-up studies with an integral held 
spectrograph. 

• A small group of galaxies is identified near the main de¬ 
hector and included in the mass model. The group of galax¬ 
ies had previously been misidentihed as a single galaxy in 
seeing limited images. 

• The redshift of the main dehector is measured to be 
z = 0.745 from stellar absorption features, revising the pre¬ 
vious estimate pf z = 0.748 based on an intervening Mg II 
absorber. 

• The stellar velocity dispersion of the main dehector is 
measured to be 290 =b 30 kms -1 . 

• Modeling the surface brightness distribution of the 
lensed source requires an accurate model of the PSF. As 
no stellar image is available in the narrow held of view of 
NIRC2 we use an iterative procedure to reconstruct the PSF 
from the lensed quasar images themselves. 

• We construct lens models of the system, including 
nearby perturbers in addition to the main dehector. The 
models reproduce the data for a close-to-isothermal total 
mass density prohle of the main dehector. Some of the lens 
model parameters suffer from relatively large uncertainties 
dominated by residual PSF uncertainties. 


• We combine the lensing and dynamical information to 
estimate the logarithmic slope of the total mass density pro¬ 
hle y . We hnd 7 ' = 1.93±0.09 consistent with typical values 
found for massive early-type galaxies at these redshifts. 

• We estimate that the current systematic uncertainties 
of the mass model of this system are too large to provide any 
interesting cosmographic inference. However, we show that, 
with the addition of HST images and spatially resolved stel¬ 
lar velocity dispersion, the system would become very useful. 
As an illustration we combine the simulated lensing likeli¬ 
hood with the Planck prior, and show that the uncertainty 
on w, curvature and Ho is signihcantly reduced. 
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